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An  Analysis  of  Physical  Optics  Models 
for  Rough  Surface  Scattering 


I.  INTRODUCTION 

There  are  a  large  number  of  approaches  to  the  calculation  of  the  electromag- 

1  — R 

netic  scattering  from  rough  surfaces.  The  importance  of  these  techniques 
rests  in  their  relevance  to  the  determination  of  the  effects  of  the  real  environment 
on  radar  and  communications  systems.  The  analyses  contain  certain  assumptions 
about  the  nature  of  the  rough  surface  in  relation  to  the  electromagnetic  phenomena. 


(Received  for  Publication  18  September  1984) 
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Particular  emphasis  has  been  placed  on  the  characterization  of  the  surface  in 

terms  of  the  statistical  distribution  of  the  heights,  their  degree  of  correlation 

2 

(T  is  the  surface  correlation  length),  the  variance  of  heights  a  ,  and  a  complex 

7  8 

dielectric  constant  representing  a  particular  terrain  type.  ’  These  features  are 

then  related  to  a  normalized  radar  cross  section  for  the  terrain  through  an  elec- 
9  10 

tromagnetic  analysis.  '  One  of  the  most  frequently  applied  approaches  is  that 
of  physical  optics. 

In  this  report,  we  will  discuss  several  topics  related  to  the  use  of  a  physical 
optics  scattering  model.  These  are  the  limitations  and  constraints  that  are  in¬ 
herent  in  the  application  of  such  a  model,  the  range  of  validity  of  the  analyses, 
how  different  representations  of  normalized  cross  section  o°  have  to  be  introduced 

for  different  parameter  ranges,  and  the  effects  of  these  different  representations 

g 

of  the  normalized  cross  section  on  the  calculation  of  the  diffuse  scattered  power. 

2.  THE  CONDITIONS  FOR  PHYSICAL  OPTICS 

Over  the  years,  the  arguments  to  support  the  validity  of  physical  optics  have 

centered  around  the  use  of  the  Kirchhoff  integral  representation  for  the  scattered 

EM  wave,  where  the  boundary  conditions  on  the  surface  have  been  satisfied 

through  the  use  of  the  Fresnel  plane  wave  reflection  coefficients.  The  boundary 

conditions  are  met  by  specifying  the  total  field  on  the  rough  surface  as  the  sum  of 

an  incident  field  and  a  scattered  field.  The  scattered  field  is  expressed  as  the 

product  of  the  incident  field  and  the  Fresnel  plane  wave  reflection  coefficient.  The 

use  of  the  Fresnel  plane  wave  reflection  cotfficient  is  justified  if  the  local  radii  of 

curvature  on  the  rough  surface  are  large  compared  to  a  wavelength  X)>  We 

wish  to  relate  this  constraint  to  more  readily  handled  surface  parameters. 

In  Ulaby  et  al^^  there  are  two  conditions  given  for  the  validity  of  the  Kirchhoff 

2 

approximation:  kT>  6  and  T  /a>  2.76  X,  where  k  =  27r/X,  X  =  EM  wavelength, 
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T  s  surface  correlation  length,  and  o  =  standard  deviation  in  surface  height.  In 
the  present  report,  it  will  be  shown  that  a  single  condition  (T»  \)  is  sufficient 
for  the  validity  of  physical  optics.  It  will  be  shown  that  for  a  rough  surface  with 
a  Gaussian  correlation  function,  T»  X  implies  that  the  radius  of  curvature 

must  be  large  compared  to  a  wavelength  (R„»  X).  Thus,  the  second  condition 

11  ^ 
given  by  Ulaby  et  al  is  redundant. 

2.1  Magnetic  Field  Integral  Equation 

To  clarify  some  of  the  implications  in  the  use  of  physical  optics  concepts  we 

4 

follow  the  treatment  of  Brown  who  starts  with  the  magnetic  field  integral  equation 
for  the  current  Jg  induced  on  a  perfectly  conducting  rough  surface  by  an  incident 
magnetic  field  H*-. 

Jjr)  =  2nXH\T)  + ^  f  n(r)X[J(r)Xv^  g(\r  -T\)]dS  .  (1) 

S  s  SOO  DO 

where  n  is  the  unit  normal  to  the  rough  surface,  S  is  described  by  z  =  |(x  .y  ), 

°  o  o  o  o 

and  g(|r  -  r^l)  is  proportional  to  the  free  space  Green's  function: 
g(lr  -  r^  I)  =  exp  (ikir  -  r^l )/  Ir  -  r"^! 


where  r^  is  the  vector  from  the  origin  to  a  point  on  the  rough  surface: 


Also,  ?(x,y)  is  the  surface  height  (a  random  variable), 


r„=x  u  +  VU  +ZU 
o  ox  *^o  y  o  z 

and  S  =  9?/9x  and  ?  =  8f/8y  are  the  x  and  y  components  of  the  surface  slope. 

^  12  ^ 

Brown  has  indicated  steps  that  transform  this  expression  into  a  new  form 
that  can  be  analyzed  to  arrive  at  a  sufficient  condition  for  the  validity  of  the  physi¬ 
cal  optics  model.  The  details  of  this  expansion  are  presented  in  the  appendix.  As 
a  result  of  those  operations  we  arrive  at  a  new  integral  equation  describing  the 
current: 


"/■/  jq<'3l''^2"  -^n’lPl<‘^l  -  ^  -  h--  -  K-  0n> 

n-H3o  T 


*4^ 


=  2H  |c^  -  i  c'l 

"  '  z  X  ok. 


q  § _ 

2x  y 


k.  ,  k.,, . .  .  k  ) 
iz'  2  n 


(2) 


12.  Brown,  G.S.  (1983)  Private  communication. 


» 


Here,  H  «  amplitude  of  the  incident  magnetic  field,  c'^,  c*^,  and  c*^  are  polariza- 
o  X  y  z 

tion  factors. 


n 


and  p^  is  the  stochastic  Fourier  transform  of  the  single  point  probability  density 
function  (pdf).  The  current  is  represented  as 


Here  k,,  k„,  .  . .  k  and  (3,,  are  stochastic  Fourier  transform  variables. 

Also,  kij.  is  the  transverse  part  of  the  incident  wave  vector  and  r^  lies  in  the 
x-y  plane.  The  details  of  these  equations  are  presented  in  the  appendix.  The 
particular  term  of  interest  in  this  expression  is  (k^,(3j,P2<  kg,  .  .  .  ,  k^ ). 


Tg  (pj,kj.P2.i^«  •  •  • 


(2 


t?//; 


3g(Ar^,  3  )  a 

a  ^x 


ag(A?,6  )  a  _  ^ 


X  Pg  (kj  -(3^,  Pq-01-  kg,. .  .kj^)exp(ik.^  '  Ar^)dAr^dj3^ 


(3) 


Here,  g  and  g^  are  the  Fourier  transforms  of  the  Green's  function  and  its  deriva¬ 
tive  respectively  and  pg  (...)  is  the  stochastic  Fourier  transform  of  the  two  point 
pdf.  The  significance  of  is  that  under  the  condition  that  this  term  vanishes 
(r^  —  0)  the  integral  equation  then  yields  the  physical  optics  result  (for  a  perfect 
conductor)  for  the  current  density: 

J  =  2  n  >  H^ 
s 


Thus  we  can  determine  a  sufficient  set  of  conditions  for  physical  optics  in 
this  instance  by  determining  a  set  of  conditions  that  causes  Tg  —  0.  The  object 
is  to  relate  dimensions  of  the  surface  to  the  wavelength  of  the  incident  field  and 
in  this  context  we  will  show  that  the  condition  T/X  -•  «  is  a  sufficient  condition  for 
Fg  0  and  hence  for  physical  optics  conditions  to  prevail. 

Since  F^  is  such  a  complicated  expression  we  simplify  the  analysis  by  exam¬ 
ining  the  type  of  behavior  expected  for  various  elements  of  the  Fg  integrals.  The 
first  term  that  we  wish  to  consider  is  pg  (.  .  . ),  the  two  point  stochastic  Fourier 
transform.  Brown^  has  given  expressions  for  PgiC ,  f  and 


4 


L'  V  iF' 


k2.  ^2’  •  •  • '  ^n’  these  relations  in  assessing  the 

behavior  of  the  multivariate  distribution  under  the  transform  operation. 

First,  we  define  t'  =  t/^  =  +  (^y)^/^  and  T'  =  T/X-.  We  restrict  our 

analysis  to  the  case  where  the  surface  correlation  function,  R(t')  is  Gaussian 

R(t')  =  CT^exp  [-  (t')^/(T')^  ] 

In  addition,  the  two  point  probability  density 

P2(C.?o'  '^0^0*  •  •  •  •  ^0^0*  ^^t' 

is  Gaussian  with  zero-mean  variates. 

Next,  we  consider  the  quantity 

<f  1^)  =  R(t')  =  exp  (  -  (t')^/(T')^  ] 

Also. 


,p  llv  -  il.'  9R<t')  . 

^  °ax  ^  '  ax  ar' 

In  the  limit  when  T'  —  oo,  (f 

'  o 

Similarly,  we  have  that  (  f  • 


R(t') 


|i)  ^  0  and  hence.  (1^ -VC)  0. 

T'.*oo 

V  C  0.  We  continue  this  evaluation  for  higher 

o’ o  aim 

T'-oo 


order  derivatives. 


(I  s!?)  .iL-ri:^  (HI) 

'  °  ax^  ax'  [  T '  \(T')2/ 


=  .  2R(t')  +  4(x'  -  x^)^R(t') 
(T')2  (tV 


This  expression  —  0  as  T'  -»  «  so  that  (  ^  ®  ^^^o^o^  lim"* 

T'-co  T'*«) 

In  similar  fashion  we  can  show  that  the  general  cases  are 

<«oV"5>  lim-  <«<^o>  Urn-  0  (  V*"!  0. 

T'-co  T'-oo  T'-oo 

In  Brown's  formalism,  the  expression  for  P2  when 

P2(E.f^.Vf,  . V"f. 

is  Gaussian  is 


P2  ^o  ■  ^2'  ■  ^2  ■  ■  ■  ■'>"  I"  ^2 


5 


where 


^2  = 


-3, 


and  C2  is  the  covariance  matrix 


C 


2 


(5^  (55o> 

<^o) 

(Vf-5N  <v5-?„) 


(  ?  •  V5  ) 

<  5o  •  V5  ) 
<(V5  )2) 

(Vo  - 


This  form  of  the  covariance  matrix  obtained  by  Brown  is  a  general  one.  For  our 
case  with  a  Gaussian  form  for  the  correlation,  we  obtain  the  simpler  form  for  the 
matrix 


■a2 

R 

0 

0 

0 

R 

a2 

0 

0 

0 

0 

0 

0 

0 

0 

^  2  ‘ 

0 

0 

0 

0 

0 

• 

• 

• 

• 

by  substituting  our  previously  derived  terms  for  the  derivatives  as  T'  -*  00  . 
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such  that  the  Rayleigh  parameter,  2  >  3.  For  smaller  values,  the  Taylor  series 
representation  for  o°  is  accurate  and  should  be  used  for  calculating  the  bistatic 
clutter  power. 

If  the  Rayleigh  parameter  dependence  of  the  cross  section  is  taken  into  ac¬ 
count  and  surface  shadow  regions  are  included,  then  these  results  indicate  that 
physical  optics  principles  can  be  applied  to  analyze  electromagnetic  scattering 
from  surfaces  with  an  extensive  range  of  roughness  conditions. 


and  surface  parameters  are  such  that  the  Rayleigh  parameter  £  ^3.0.  For 
smaller  Rayleigh  parameters,  the  Taylor  series  representation  for  a*  is  accurate 
and  should  be  used  for  calculations  of  clutter  power. 

i  CONCLUSIONS 

In  this  report  four  major  topics  have  been  studied.  The  findings  for  each  are 
summarized  here.  The  first  topic  is  the  investigation  of  the  conditions  under 
which  a  physical  optics  model  of  the  scattering  from  a  rough  surface  is  a  valid 
representation.  The  investigation  proceeded  from  an  analysis  of  the  magnetic 
field  integral  equation  for  the  current  density  on  the  rough  surface.  For  a  surface 
with  a  Gaussian  surface  height  distribution  and  a  Gaussian  surface  correlation 
function  we  have  demonstrated  that  only  one  sufficient  condition  is  needed  in  order 
that  physical  optics  be  valid  for  calculating  EM  wave  scattering  from  the  rough 
surface.  This  condition  is  that  the  correlation  length  T  be  much  greater  than  a 
wavelength. 

The  second  consideration  of  the  report  is  to  relate  that  constraint  to  the  al¬ 
ternative  condition  that  the  radius  of  curvature  of  the  surface  irregularities  be 
large  compared  to  the  wavelength.  The  analysis  involved  the  same  restrictions  of 
Gaussian  surface  heights  and  surface  correlation  function.  For  our  condition 
(T/X—  00  )  it  was  shown  that  the  consequence  of  these  assumptions  is  that  the  sur¬ 
face  slopes  and  slope  derivatives  are  statistically  independent.  The  conclusion 
that  follows  from  these  arguments  is  that  T  »  X  is  always  a  sufficient  condition 
for  X  and  for  most  cases  it  is  also  necessary. 

In  the  third  topic  we  have  addressed  the  question  of  the  form  of  the  scattering 
cross  section  that  is  present  in  the  physical  optics  representation  of  the  scattering. 
In  the  report  we  have  shown  that,  despite  the  constraint  on  surface  correlation,  the 
surface  can  be  described  by  a  wide  range  of  values  of  the  Rayleigh  roughness 
parameter.  Typically  in  radar  engineering,  the  asymptotic  form  for  the  cross 
section  is  used  in  scattering  calculations.  Here  we  have  shown  that  different  rep¬ 
resentations  for  the  normalized  cross  section  a"  must  be  used,  depending  on  the 
magnitude  of  the  Rayleigh  roughness  parameter.  Explicit  representations  for  o' 
have  been  given  for  small,  intermediate,  and  large  Rayleigh  roughness  parameter 
when  the  surface  is  described  by  either  a  bivariate  Gaussian  distribution  or  a  bi¬ 
variate  exponential  distribution  in  height. 

Finally,  we  examined  the  effect  of  using  the  different  representations  for  o° 
in  determining  the  diffuse  power  scattered  in  the  forward  direction  by  a  bistatic 
radar  system.  The  conclusion  is  that  the  conventional  asymptotic  form  for  a° 
gives  accurate  results  (±  1  dB)  only  when  the  system  and  surface  parameters  are 
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Table  3.  Diffuse  Scattered  Power  for  Different  a*  Representations 


Antenna 

Separation 

a2 

Strip 

Distance 

Diffuse  Scattering 

Power 

92  km 

im^ 

10  m 

531. 

PDIFF  =7.52.10-13w 
PDIFFI  =7.57.10-13w 
PDIFFX  =  2.  1.  10-113W 

1  km 

6.67 

PDIFF  =3.88.  10-12W 
PDIFFI  =  4.  16.  10-12W 
PDIFFX  =  4,  16.  10-12W 

4  km 

0.795 

PDIFF  =1.09.  10-12W 
PDIFFI  =6,75.  10-13W 
PDIFFX  =  6.75.  10-13W 

46  km 

0.4 

PDIFF  =7.28.  lO-lSw 
PDIFFI  =4.41.10-15W 
PDIFFX  =  4.41.  lO-lSW 

10  m^ 

46  km 

4.3 

PDIFF  =3.65.  10-14W 
PDIFFI  =3.52.10-14w 
PDIFFX  =  3.  52.  10-14W 

0.  5  m^ 

46  km 

0.2 

PDIFF  =5.  53.  10-15W 
PDIFFI  =2.27.10-15W 
PDIFFX  =  2. 27.  10-15W 

0.  1  m^ 

46  km 

0.  04 

PDIFF  =2.67.10-16W 
PDIFFI  =4.67.10-16w 
PDIFFX  =  4.67.  10-16W 

10  m^ 

10  m 

5308. 

PDIFF  =6.75.  10-12W 
PDIFFI  =6.75.  10-12W 
PDIFFX  =2.08.  10-112W 

9  km 

10  m^ 

10  m 

6663. 

PDIFF  =3.04.10-9W 
PDIFFI  =3.04.10-9W 
PDIFFX  =  2.08.  10-1 13W 

_ 1 

1  km 

328. 

PDIFF  =3.62.10-13W 
PDIFFI  =3.63.10-13W 
PDIFFX  =  2.08.  10-113W 

contribution  for  two  strips  have  been  evaluated,  one  close  to  the  transmitter  and 

2 

the  second  at  a  distance  from  it.  These  comparisons  show  the  effect  of  2  on  the 
three  types  of  power  calculations  as  a  number  of  parameters  are  allowed  to  vary 
widely. 

In  calculations  of  radar  clutter  contributions  from  Gaussian  surfaces,  the 
asymptotic  form  of  a°  is  conventionally  used.  Examination  of  the  results  of  this 
table  indicates  that  the  asymptotic  form  (physical  optics  plus  steepest  descent 
evaluation  of  integral)  for  a°  will  give  accurate  results  (±  1  dB)  only  if  the  system 


4.  EFFECT  OF  a°  REPRESENTATION  ON  DIFFUSE  SCATTERED  POWER 


In  this  section,  actual  quantitative  comparisons  will  be  made  of  the  diffuse 

power  reaching  a  receiver  (see  Figure  1),  where  the  different  representations  for 

CT°  have  been  used,  as  discussed  in  Section  3.  Because  the  computation  time  was 

prohibitive  when  the  total  diffuse  power  was  calculated  by  integrating  over  the 

entire  glistening  surface  (especially  for  the  integral  representations  for  a°),  we 

consider  only  the  diffuse  power  reaching  the  receiver  from  selected  strips  of 

ground  20  m  wide  and  parallel  to  the  y-axis. 

The  effect  of  integrating  across  the  surface  can  be  seen  in  the  results  shown 
2  2 

in  Table  3.  For  a  =  1  m  and  an  antenna  separation  of  92  km,  we  show  the  in¬ 
dividual  contributions  from  strips  that  represent  a  wide  variation  in  position  along 
the  distance  axis.  We  are  concerned  with  the  effect  of  the  Rayleigh  parameter  2 
on  the  various  representations.  There  are  two  main  factors  constituting  that 
parameter,  a  and  scattering  angle.  In  the  table,  ^  results  are  obtained,  where 
both  of  these  factors  are  allowed  to  vary  separately.  The  results  presented  here 
are  for  a  receiver  height  of  2  500  m,  transmitter  height  of  100  m,  transmitted 
power  of  350  W  and  vertical  polarization.  The  surface  heights  are  taken  to  have  a 
Gaussian  distribution  and  a  surface  correlation  length  T  =  5  m.  For  all  cases,  we 
show  the  three  different  diffuse  power  levels  calculated  using  the  three  represen¬ 
tations  of  O’";  PDIFF  corresponds  to  the  asymptotic  form;  PDIFFI  is  the  numerical 
evaluation  of  the  integral;  and  PDIFFX  is  the  result  for  the  power  series  repre¬ 
sentation. 

There  are  three  basic  comparisons  in  the  table.  First,  for  a  separation  of 
2  2  2 

92  km  and  a  =  1  m  the  effect  of  2  is  shown  as  the  distance  of  the  scattering 

strip  (scattering  angle)  is  varied.  Next,  for  the  same  separation  (92  km)  and  a 

2  2 
fixed  strip  position  (46  km  from  the  transmitter)  2  is  changed  by  varying  o 

2  2  2 

(0.  1  m  s  a  10  m).  Finally,  the  table  shows  a  comparison  of  results  for 
2  2 

a  =  10  m  at  two  antenna  separations  (9  and  92  km).  At  both  separations  the 


ROUGH  terrain 


Figure  1.  Rough  Surface  Scattering  for  a  Bistatic  Geometry 


Table  2.  J°  Values  Calculated  Using  Different 
Representations 


J*  for  Bivariate  Gaussian 

z2 

POWER  EXP. 

ASYMPT. 

INTEGRAL 

96 

1.  OE  -  13 

7807 

20 

7323 

8018 

10 

8300 

7753 

8298 

1 

3850 

77  53 

9931 

J°  for  Bivariate  Exponential 

^2 

POWER  EXP 

ASYMPT. 

INTEGRAL 

96 

9881 

9850 

10 

155779 

9881 

9985 

1 

7753 

9881 

3735 

0.  1 

91 

9881 

85 

The  results  of  this  table  may  be  summarized  as  follows: 

(  i)  Bivariate  Caussian  Surface: 

2 

(a)  Tor  large  Rayleigh  parameters  >  20),  both  the  asymptotic  and 

integral  representations  are  accurate  and  agree  to  within  a  few  percent. 

2 

(b)  For  intermediate  and  large  Rayleigh  parameters  (2  >  10),  the  in¬ 

tegral  representation  is  accurate  (for  small  ^  values  the  numerical  integration 
may  have  oscillation  problems). 

2 

(c)  For  small  Rayleigh  parameters  (S  «  10),  only  the  power  series 
expansion  representation  is  accurate. 

(ii)  Bivariate  Kxponential  Surface: 

2 

(a)  For  large  Rayleigh  parameters  2  20),  both  the  asymptotic  and 
integral  representations  are  accurate  and  agree  to  within  a  few  percent. 

(b)  For  all  values  of  the  Rayleigh  parameter,  the  integral  representation 
IS  accurate 

2 

(c)  For  small  values  of  the  Rayleigh  parameter  (C  <  1),  the  power 
series  expansion  representation  is  accurate  and  agrees  with  the  integral  repre¬ 
sentation  to  within  a  few  percent. 


17 

Brown  has  identified  Jp  as  the  incoherent  power  scattered  diffusely,  and 
Jg  is  the  incoherent  power  scattered  in  the  specular  direction  (v^  ” 
four-fold  integral  given  in  Eq.  (9a)  can  be  reduced  to  a  single  integral  if  T  «X,  Y, 
just  as  was  true  for  the  Gaussian  case.  The  result  is 


■>0 


(10) 


2  D 

For  large  and  small  values  of  ,  J  ==  0  so  J°  ==  J  .  For  large  Rayleigh  para- 
2  ^ 

meters  (E  >  20),  a  steepest  descent  evaluation  of  the  integral  in  Eq.  (10)  yields 
the  asymptotic  expression  for  J^: 


(11) 


2 

where  =  -  (2  7r/X)(cos  0^^  +  cos  0^).  For  intermediate  values  of  2  ,  the  integral 

has  to  be  evaluated,  and  in  general,  J  ^0  so  J*  =  J^.  +  J  .  For  small  Rayleigh 

2  ^  ^ 
parameters  (S  <  1),  the  following  expression  [from  a  Taylor  expansion  of  the 

integral  in  Eq.  (10)]  can  be  used; 


J 


D 


E^T^  exp 


2  ^2 
-  V  T 

_ 5:1: _ 


(12) 


.3.3  CompariKiiu 

To  demonstrate  quantitatively  the  accuracy  of  these  various  representations 

2 

for  J*  as  a  function  of  the  Rayleigh  parameter  squared  (2  ),  Table  2  was  con- 

2 

structed.  This  table  shows  J*  values  as  a  function  of  ^  using  the  different  rep¬ 
resentations.  Here,  POWER  EXP.  refers  to  the  Taylor  series  expansion  repre¬ 
sentations  [Eqs.  (8)  and  (12)|,  ASYMPT.  refers  to  the  asymptotic  representations 

[Eqs.  (7)  and  (11)|,  and  INTEGRAL  refers  to  the  integral  representations  [Eqs. 

2  2 

(6)  and  (10)|.  In  these  tables,  the  variance  in  surface  height  a  =  10  m  ,  the  cor¬ 
relation  length  T  =  1.^5  m,  and  the  azimuthal  scattering  angle  *  0.  5*.  It  should 
be  noted  that,  since  we  are  concerned  with  differences  between  solutions,  the 
values  shown  for  the  intermediate  exponential  solutions  are  actually  values  for  Jj^ 
since  Jg  would  be  nearly  zero  for  all  methods. 


17.  Brown,  G.S.  (1982)  Scattering  From  a  Class  of  Randomly  Rough  Surfaces, 
Radio  Science  17:1274-1280. 
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Next,  we  consider  the  evaluation  of  this  expression  for  different  2  regimes.  If 

2»  1,  a  steepest  descent  evaluation  of  the  integral  in  Eq.  (6)  reduces  it  to  the 

1 6 

asymptotic  expression  familiar  to  radar  engineers  (Barton  and  Ward  ): 


Ia  2-2^ 

r-v2  T^l 

/4s-  T 

jexp 

__xv 

J 

1 

M 

1 _ 

(7) 


In  the  intermediate  case,  the  integral  has  to  be  evaluated  explicitly.  When  the 

2 

Rayleigh  parameter  is  small  (2  <  10),  the  following  series  representation  can 

be  used 


J°  = 


f4,Vl 

v2  " 

22m 

exp 

r  2  ,„2‘i 

-V  T 
xy 

\  i 

(m.  m! ) 

4m 

m*l 

- 

(8) 


3.2  Bivariate  Exponential  Solution 

The  preceding  results  are  all  for  Gaussian  distributed  surface  heights.  For 
the  case  where  the  surface  heights  are  described  by  a  bivariate  exponential  the 
four-fold  integral  appearing  in  Eq.  (5)  may  be  written  as  follows: 

J°  =  J  +  J,,  ,  (9) 

o  s 

where 

X  X  y  y 

/dYg  exp[iv^(Xj  -X2)+iVy(y,  -y^)) 
^  “X  “X  “y  “y 

[Xg  -  \2(t  -oc)]  (9a) 

and 

Js  =l-^j(4XY)  sinc^  (v^X)  sinc^  (v^Y)|  (1  +|- 22)-3'2  .  (j  +.1  v;2)'2]  (pb) 

Since  this  is  for  an  exponential  distribution,  the  bivariate  characteristic- 
function  for  uncorrelated  variates  used  in  Eq.  (9a)  and  Eq.  (9b)  has  the  form 

X  2(t  -oo)  =  [  1  +  (2/3)22i-2/2 

16.  Barton,  D.K.,  and  Ward,  H.R.  (1969)  Handbook  of  Radar  Measurement, 
Prentice -Hall,  Inc.  Englewood  Cliffs,  New  Jersey. 


XYX" 


A  *  y  J- 

J  '^1  /dY.  y* dYg  exp[iv^(X^  -Xg)  +  iVy(Yj  -  Yg)] 

J  -X  -X  -y  -y 


X[X2(2:.T)  -  Xi(2)Xj(£)l 


Here,  are  the  matrix  elements  for  linear  polarization  states  (see  Ruck  et  al'^), 
the  rou^  surface  has  an  area  4XY,  and 


T  =  Surface  Correlation  Length 

v„  =  (2  ff/X)[  sin  e.  -  sin  0  cos  (^  ] 

X  I  s  s 

V  =  -  (2:r/x)  sin  0„  sin  </)„ 

J  o  a 

X  j  =  univariate  characteristic  function  of  the  surface  height 
distribution  function, 

X  2  =  bivariate  characteristic  function, 

2  =  Rayleigh  roughness  parameter, 

=  (2  Tr<y/X)  (cos  0.  +  cos  0„  ) 

1  s 

The  forms  assumed  by  x  ,  and  Xn  for  Gaussian  surface  height  distributions 

^  ^  9 

and  for  exponential  surface  height  distributions  have  been  derived  previously. 

For  Gaussian  surfaces,  Xg^’’’  )  =  sxp  (  -  2  )  and  Xg  -•  » )  =  x  j,i,X  j-  On  the 

other  hand,  for  a  surface  described  by  a  bivariate  exponential  distribution  function, 

2  •3/2 

X 2  (t  —oo)  =  { 1  +  (2/3)  2  j  '  and  Xg  (t  — oo)^  X  j;»X  j.  The  only  condition  needed 
to  reduce  the  four-fold  integral  to  a  single  integration  is  T  «X,  Y.  This  criterion 
does  not  depend  on  the  Rayleigh  parameter  but  there  are  differences  in  the  result¬ 
ing  forms  for  the  two  surface  distributions. 


3.1  Gaussian  Solutions 

For  the  Gaussian  case,  we  have; 


2 


(6) 


where 


hence  in  that  range  T  »  \  »  X  and  R  »  X  -►  T  »  X.  Thus,  T  »  X  is 

both  necessary  and  sufficient  for  R^  »  X  when  we  have  intermediate  o/T  values. 

To  summarize  the  above  results,  we  can  say  that,  under  our  assumptions  of 
independent  slope  and  slope  derivatives,  for  a  Gaussian  surface  and  Gaussian  cor¬ 
relation  function.  T  »  X  is  always  a  sufficient  condition  for  R^  »  X,  and  for 
intermediate  and  large  slope  conditions,  T  »  X  is  both  necessary  and  sufficient 
for  a  physical  optics  solution  to  apply. 

2.3  Shadowing 

There  is  a  final  aspect  to  the  relationship  between  roughness  and  physical 
optics.  It  is  well  known  that  as  the  average  surface  slope  becomes  large  (ct/Ts  1), 
shadowing  becomes  more  and  more  important.  Then,  the  physical  optics  current 
density  J  =  2(n  x  H  does  not  truly  hold  everywhere  on  the  surface.  For  large 
slopes,  the  correct  current  is  obtained  by  multiplying  the  a°  resulting  from  physi¬ 
cal  optics  by  a  shadowing  function  S,  which  describes  how  much  of  the  surface  is 
unlit,  that  is,  where  J  —  0  in  the  shadow  regions.  At  present,  this  shadowing 
function  has  been  derived  rigorously  only  in  the  high  frequency  geometrical  optics 
limit  (  X  -•  0). 

3.  NORMALIZED  SCATTERING  CROSS  SECTION 

7-10 

Previous  reports  and  papers  discussed  how  a  general  expression  for  the 
normalized  cross  section  a”  may  be  derived  from  the  Kirchhoff  integral  expres¬ 
sion  for  the  waves  scattered  from  a  rough  surface.  In  the  present  report,  a  brief 
outline  will  be  given  for  the  determination  of  <7°  for  different  regimes  of  the  Ray¬ 
leigh  parameter.  The  actual  evaluation  of  the  four-fold  integral  depends  upon  the 
surface  height  distribution  function.  The  general  expression  is 


14.  Sancer,  M.I.  ( 1969)  Shadow -corrected  electromagnetic  scattering  from  a 

randomly  rough  surface,  IEEE  Trans.  Antennas  Propag,  AP-17 :577  -585. 

15.  Brown,  G.S.  (1980)  Shadowing  by  Non -Gaussian  random  surfaces.  Proceed¬ 

ings  of  the  Second  Workshop  on  Terrain  and  Sea  Scatter.  George  Washing¬ 
ton  University,  Washington,  D.  C. 


y  ^  '  J*.  ,  i”.  .r,  r. 


Thus.  T»  \  is  both  necessary  and  sufficient  for  X  when  ff/T»  1. 

In  the  intermediate  range  of  X  the  limiting  case  forms  do  not  apply  and  the 

function  U(l/2,  0.  X)  must  be  evaluated  explicitly.  The  procedures  for  evaluation 

13 

are  outlined  in  Abramowitz.  To  see  how  varies  in  the  intermediate  X  range 

we  consider  several  cases  as  10.0  >X  >0.  1.  Over  this  range,  1.  5  2(a/T)  >0.  16 

and  the  evaluation  of  U(l/2,  0.  X)  leads  to  the  result  that  0.  5T  £  R  £  2.  5T.  In 

c 

the  earlier  report  we  evaluated  R^  under  the  assumption  that 

WW/ ■  <vr®)  ' 

If  we  compare  results  of  that  early  approach  to  those  of  the  present  case  for  inter¬ 
mediate  slope  conditions  we  see  that  for-a/T<  1  the  results  are  in  reasonable 
agreement  while  they  diverge  for  o/T>  1.  These  results  are  summarized  in 
Table  1. 


Table  1.  Results  for  in  the  Intermediate  Range  of  o/T  Values 


X 

a/T 

(|K|> 

Rc 

Previous  Rj, 

10 

0.  16 

0.4/T 

2.  5T 

1.  9T 

3 

0.3 

0.  5/T 

2.  IT 

1.2T 

1.  6 

0.4 

0.7/T 

1.5T 

1.  IT 

1 

0.5 

0.9/T 

1.  IT 

1.  IT 

0.  5 

0.7 

1.  1/T 

0.  9T 

1.  15T 

0.4 

0.  8 

1.  3/T 

0.  8T 

1.24T 

0.  1 

1.  5 

1.9/T 

0.  5T 

2.7T 

It  should  be  noted  that  in  a  strict  sense  the  regions  of  X  and  a/T  are  not 
equivalent.  For  instance  a/T  =  10  —  X  =  0.0025.  which  is  well  into  the  small  X 
approximation  and  similarly  a/T  =  0.  1  -*  X  =  25,  which  is  in  the  large  X  solution 
region.  Equivalently,  for  0. 1  s  X  s  10  we  have  1.  5  >  a/T  >  0.  16.  Thus,  for 
completeness  we  should  examine  the  small  X  solution  for  1.  5  s  a/T  ^  10  to  com¬ 
plete  the  examination  of  how  R^  behaves  for  intermediate  a/T.  For  that  range 
R  =  0.6T.  Similarly,  for  0.  16  >  a/T  >  0.  1  we  have  2.3  T  £  R  £  3.6T.  Thus 
over  the  entire  range  of  intermediate  a/T  values  R^  is  of  the  same  order  as  T  and 


13.  Abramowitz,  M. ,  and  Stegun,  I.  A.  (Ed. )  (1964)  Handbook  of  Mathematical 
Functions  With  Formulas.  Graphs  and  Mathematical  Tables.  NBS  Applied 
Mathematics  Series  No.  55. 
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After  some  manipulation  we  obtain 


where  U(a.b.X)  is  the  confluent  hypergeometric  function  of  the  second  kind  with 
X  s  0.5 

We  now  want  to  examine  the  relationships  for  various  slope  regimes.  Recall 
that  (z'^)  =  2o^/T^  and  (z"^)  =  12a^/T^.  Then,  for  small  slopes  (o/T  «  1)  we 
have  (2  (z'^)  ^  —  oo  and 


{|K|)  = 


(^?)  [■P'\ 


-1  / _ 1_\  -1/2 

2(z*2) 


)  s  2.76a 


for  zeroth  order  and 
2.76a 


<lKl) 


1  -  3 

a] 

2' 

^Tj 

for  first  order. 

Thus,  we  have  a  zeroth  order  solution  for  R  s  |  ) 

•=  \<|K|>/ 


c  2.76a 


which  is  equivalent  to  the  result  obtained  in  Ulaby  et  al^^  for  small  slopes.  Then, 
since  a/T  «  1  we  have  R^^  0.36T  and  T»  X  (a  sufficient  condition). 

For  large  slope  conditions  (a/T»  1),  /  — ^  0  and 

l2(^'  >j 


so 

R(.  s  0.  6T 
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Since  they  are  considering  only  small  slopes,  the  curvature  is  considered  to  be  a 
function  of  just  the  slope  derivatives.  For  our  applications,  though,  we  are  con¬ 
cerned  with  behavior  for  all  surface  slope  magnitudes  and  therefore  we  cannot  use 
their  simplified  expression  for  the  curvature.  As  a  result,  the  analysis  is  more 
complicated. 

If  we  make  use  of  the  previously  cited  theorem,  we  can  express  the  proba¬ 
bility  density  function  for  the  surface  slopes  as: 


and  for  the  slope  derivatives 


/  1  \ 

z"2 

_  2<z"2) 

In  the  previous  section  we  showed  that,  if  the  surface  correlation  length  T'  —  »  , 
then  the  slopes  and  slope  derivatives  are  statistically  independent  with  (z'z")  -*•  0. 
Hence,  under  our  assumptions  we  can  write  the  joint  probability  density  p(z',  z") 
as 


p(z',z")  = 


V 

-J  exp 

[li 

,.2 


.  "2 


2r(z'^)  ) 

Next,  we  rewrite  the  curvature  magnitude  as 

IkI  =  |z"[  1  +  =  Iz"l  1  [  1  + 

Then  for  the  expected  value  we  have 


..2^ 


,ti2. 


nKl' 


lz"l  1  [  1  +  z'^]’^/2j  p(z',z")dz'dz" 


Integrating  over  z"  results  in 

^iKh = 


•j  ^  [  1  +  z'^]  exp  [  -  z'^/(2(  z'^N )]  dz' 
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R  is  the  average  radius  of  curvature.  This  is  the  most  generally  used  form  for 
the  autocorrelation  since  use  of  an  exponential  form  leads  to  discontinuities  at 
small  distances. 

g 

In  a  •previous  report  by  Papa,  Lepnon,  and  Taylor,  a  derivation  was  given 
for  the  mathematical  representation  of  the  average  radius  of  curvature  R^.  The 
rough  surface  is  described  by  the  equation  z  =  ?  (x,y).  The  average  slope  is  given 

by 

.  r  /a,  1^/2 


where  (  )  denotes  expectation  value.  The  average  radius  of  curvature  R^  of  a 
curve  is  given  by 


^C  “  <|k|) 


[(1  -H  (Z')^]  ^^^1 


where  K  is  the  curvature.  The  average  of  {z'}  is  given  by 


.  /«  !iV  -  ■  K 

\  ax  ax  /  a 


Also,  the  average  of  (z'^r  is  given  by 

.  !2|!  , 

a  T  ^  t"* 

The  results  are  consistent  with  the  theorem  from  random  processes  that  states 

that  the  distribution  of  the  derivative  of  a  normal  process  with  zero  mean  and 
2 

variance  a  is  again  normal  and  the  variance  of  the  derivative  is 

Og  =  -o  p"  (0)  . 

This  theorem  will  be  used  by  us  to  examine  the  relationship  between  the  radius  of 
curvature  and  the  correlation  length  in  physical  optics  models. 

Ulaby  et  aK*  examined  this  relation  for  the  case  where  the  slopes  are  re¬ 
stricted  to  the  condition  that  (z')  «  1.  They  obtain  the  result  that 


2.76a 
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Using  this  simplified  matrix,  we  obtain  the  result  that  in  the  limit  as  T'  oo 
^2^^!  *  ^o'  ^1  '  ^o'  ^2’  ~  ^2'  ■  ■  •  ^ 

-  exp  (- 4  {(ki  -  3^)^  +  2(3^  - 3i)(kj  - 3i) +  (3o  - 0i)^  }  1 

These  results  now  allow  us  to  examine  the  individual  terms  of  in  Eq.  (3). 
Since  pg  is  independent  of  in  the  limit  as  T*  -*  oo  the  terms 

I  dgi^r^,  3q)  \  and  /  1  vanish.  To  examine  the  final 


term,  we  reformulate  the  integration  as 


00 

1=  //y/i3-g(?)exp[-i3^?)P2  X 
-00 

dAx  dAy  d^^  df  , 


and  make  an  asymptotic  evaluation  of  the  Ax  and  Ay  integrations,  using  the 
stationary  phase  method  with  k^  =  k^T  -*  oo  .  This  leaves  a  double  integration  in 
?  and  3q-  The  3^  terms  in  the  integrand  that  results  from  the  asymptotic  evalu¬ 
ation  are  a  standard  form  so  that  the  entire  process  then  reduces  to  a  single 
integral  in  the  complex  I -plane.  Careful  analysis  of  the  appropriate  contour 
integration  shows  that  the  integral  reduces  to  zero.  The  final  result  of  this 
analysis  is  that  for  T/X  -*  w.  ^2  —  0.  This  implies  that  Eq.  (2)  reduces  to 


^  /jq^01>Pl<ki  -3i)d3i  =  2c^?j(kj  -k.^) 


with  q  =  X  or  y.  The  solution  of  Eq.  (4)  is  given  by  the  result 
jq(z)  =  2H^  c^  exp  ),  which  is  the  physical  optics  current.  Thus,  T/X  — oc 

implies  that  the  current  density  on  a  rough  surface  is  given  by  the  physical  optics 
result  [  J  =  2n  X  H  (for  a  perfectly  conducting  surface). 

2.2  Relation  Between  R^  and  T/  X 

We  have  derived  the  physical  optics  current  from  the  magnetic  field  integral 

equation  using  T/x  -»  oc  .  We  next  show  the  relation  between  this  condition  and  the 

basic  condition  for  physical  optics,  namely  that  R^»  X.  We  will  show  that  for  a 

rough  surface  with  an  autocorrelation  function  R  »  e“''^/T^,  where 

2  2  1  /2 

T  =  ((Xj  -  Xg)  +  (yj  -  yg)  ]  '  ,  the  condition  T»  X  implies  R  »  X,  where 
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Appendix  A 


An  Analysis  of  the  Magnetic  Field  Integral  Equation 


A 1 

Brown  has  outlined  how  a  general  expression  of  the  magnetic  field  integral 

equation  for  the  current  J  induced  on  a  perfectly  conducting  rough  surface  by  an 

s 

incident  magnetic  field  Hi  can  be  expanded  into  a  form  that  will  allow  us  to  analyze 
conditions  sufficient  for  a  physical  optics  model  to  be  applied  to  the  scattering 
from  a  rough  surface.  The  details  of  the  analysis  are  presented  here. 

The  current  is  given  by: 


J  (r  )  =  2n  X  hVt  )  +  /"  n(r  )  X  [  J„(r  )  x  v„g  ( I  r  -  r  I )]  dS  ,  (Al) 

o  J  oOO  OO 

where  n  is  the  unit  normal  to  rough  surface  S  described  by  z  =  §  (x  ,  y  )  and 
g(  i  r  -  r^l)  is  proportional  to  the  free  space  Green's  function: 

iTlF)  =  [  -  lyff- F]  [l  +  4  +  €yl'^/^ 

g(  1  r  -  r^  1 )  =  [  exp  (ik  1  r  -  r^  I )  1  [  I  r  -  r^  |  ]  '^ 


Al.  Brown,  G.  S.  (1982)  A  stochastic  Fourier  transform  approach  to  scattering 
from  perfectly  conducting  randomly  rough  surfaces,  IEEE  Trans.  Anten- 
nas  Propag.  AP-30  (No.  6):1135-n44. 


where  r  Is  the  vector  from  the  origin  to  a  point  on  the  rough  surface: 


'  *o"x  yo“y  *o'‘z  +  y«“. 


OX  •'o  y 


Also.  i(x,  y)  B  z(x.  y)  is  the  surface  hei^t  <a  random  variable)  and  €  > 

d€  X  aX 

and  ly  B-gA  are  the  x  and  y  components  of  the  surface  slope. 

To  examine  the  terms  in  the  magnetic  field  integral  equation  [Eq.  (Al)],  the 
double  cross  product  is  expanded,  the  surface  integration  is  converted  to  an  inte¬ 
gration  over  the  z^  B  0  plane  through  dS^g  dr^^  and  both  sides  of 

the  equation  are  multiplied  by  V  1  +  ^y 


3(r)  =  2frxfl"(?)  +  (2F)'^y’{(f}(?)  •  Vj,glT(?^) 


-  (N(r)  •  J(r^)l  V^gldr^ 


where 


lir)  = 


15  (r)  =  (7  1  +  ?2  +  52)-(-) 
X  y 


and 


dr,  B  dx  dy 
to  o  ■'o 


The  current  must  always  be  tangential  to  the  surface,  so  that 


N  (r  )  •  T(r  )  *  0 
o  o 


and  therefore. 


J-O  =  J  <rJ  +  ?..  J  (?  ) 


z  o 


X  X  o  y  y  o 
o  •'o 


(A2 


(A3 


If  Eq.  (A3)  is  substituted  into  Eq.  (A2),  the  result  will  give  two  coupled  inte¬ 
gral  equations  for  J.^(r)  and  Jy^r).  The  equations  will  become  uncoupled  if  the 
term  [N  •  J  Ir^)]  in  Eq.  (A2)  is  ignored.  This  is  equivalent  to  neglecting 
cross  polarization  contributions  to  the  current.  These  effects  could  be  included 
in  the  analysis  but  it  becomes  extremely  cumbersome. 

In  general,  the  electric  field  in  the  far  field  may  be  expressed  in  terms  of  the 
current  density  as  follows; 


E 


(A4) 


where 


r  =  xu^  +  yUy  +  Cu^ 

k=k  u+k  u+k  u 
s  sx  X  sy  y  sz  z 

dF^  =  dx  dy 
r)^  =  impedance 


and  is  the  distance  from  the  z  =  0  plane  to  the  observation  point. 

The  statistical  moments  of  the  scattered  field  (E^)  and  (  can  be  deter¬ 

mined  from  Eqs.  (A2),  (A3),  and  (A4),  where  (  )  denotes  an  ensemble  average.  It 
is  possible  to  obtain  an  equation  for  the  current  in  the  Stochastic  Fourier  Trans 
form  domain  by  multiplying  Eq.  (A2)  on  both  sides  of  the  equation  by  the  Fourier 
kernel 


00 

exp  I  -  ik  ?  -  i  ^  k^^  ■  1 

n=l 


and  then  averaging  over  all  stochastic  variables  using  the  joint  probability  density 
function  relating  all  the  h'^ights,  slopes,  and  slope  derivatives.  Performing  these 
operations  on  the  left-hand  side  of  Eq.  (A2)  yields  the  result: 


<J  (F)exp[-  ikj?  -  i  ^  •  v"f]> 


n=  1 


=  lim  (2jr)'^n /•  •  •/  J  (F,0j,02*  ••• 


Pi  (ki  -  3i,  ^2"  ^2*  *  '  ’ '  *^^2  '  '  * 


where 


n 


n 


=  E  i  . 


i=l 

r  »  xF  +  yF 


^  «X/ 

Jqis  the  Stochastic  Fourier  transform  of  Jq(r),  q  *  x  or  y  and  Pj^(. . . )  is  the 
Stochastic  Fourier  transform  of  the  joint  probability  density  function  (p  d  f).  At 
this  point  we  have  introduced  the  uncoupled  current  terms. 

If  the  incident  magnetic  field  of  the  surface  is  assumed  to  have  the  form 

H^(r)  =  exp  (+ ik.^  . 

then  the  term  2q  •  [NX  h’)  on  the  right-hand  side  of  Eq.  (A2)  may  be  written 

2q.  [N(r)xH\r)l  =  2H^lcJ+  c^C^+  cj?yl  exp(  +  iki^  .  r^  +  iku5) 

The  polarization  factors  are  determined  by  the  direction  iTj^  of  the  incident  mag¬ 
netic  field; 

Cz  =  > 

c^  =  -  q  .  X  n„  ) 

=  -  q  •  <Uy  X  njj  ) 


The  Stochastic  Fourier  transform  of  this  term  may  be  written: 


{2q  ■  (N  X  H.)  exp[  -  ik^C  -  i  k^^^  •  v"C  ]> 

n=l 

■‘^iz'^'S  •  •  ' 

(A6) 


2H  exp(+ik.. 

O  It 


^t) 


q  .  q  a 

c  -  ic 
z  X  9k, 


-  ic 


q  d 


2x 


y  9^2, 


where 

S  '  '^x  *^x  '^y^y 

As  before,  neglecting  the  coupled  term  [N  •  J(r)]  g  on  the  right-hand  side 
of  Eq.  (A2),  the  Stochastic  Fourier  Transform  of  the  integral  in  Eq.  (A2)  may  be 
written: 


( 


I 


I 


3 
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oo 

/[-  -  ‘  S  ^n+l  •5”« 


dr  V 


f-  ■  f  '  ^  ^  ■*■  If  '^y^r.'  ■  •  • 


“  ik 


11=1 


P2(?.4„.V|,V^Io-  ••• 


(A7) 


where  P^  is  the  two  point  pdf  of  the  variates.  The  £ -integration  ran  be  written  as 
a  convolution  of  the  ? -Fourier  transforms  of  the  Green's  function  derivatives  and 

P2(...). 


By  using  the  relations 

=/  g  (Ar^,  £)  exp  ( +  i3^£)  dC 


and 


«?<^^t'^o>  ■/ 


3g(Arj,?) 


31 


exp  (i3Qi)  d£ 


and  making  the  substitution  Ar^  "  ‘  ^to’  written 


< 


/[■«ylf-^-5y5f-/5f„ 


exp  -ik^?  -i3^(?^-£)-iX: 


n=  I 


drt  > 
o 


3g(Ar  3„)  3g(^r  3  )  ^  ^ 

_ 1 2.  -  c  - 1 — 2-  ^A,- 


X  3x 


-  I 


-g^(Ar^.3^) 


q  Iq  o  o  o 


u. 


K  exp 


n*l 


P2<‘^l-'^o'  V  Vo'--->'^V'^o'^^c 


(AS) 


where  P2  is  the  Fourier  transform  of  P2(.  •  •  )  with  respect  to  £. 

The  integration  over  £  may  be  expressed  as  a  convolution  with  a  shifted  argu- 
°  2 

ment  due  to  the  exponential  factor  in  Eq.  (A8).  The  ?p,  •  •  ■  integrations 

may  be  written  as  convolution  with  no  shift  in  the  argument,  and  the  integrations 
2 

over  V?,  V^.  •  •  •  may  be  written  as  simple  Fourier  transforms. 

If  these  transforms  are  performed,  Eq.  (AS)  becomes 
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I 


■  ■  ‘^ok,  ■ 

X  '^q^^'t^'^r  ^2*  ■  ■  ■  ®n^^2^^l"®o*^o*®l*''2*  '^2  '  *  ‘  ‘  ‘  ‘  **^11 


where 


@1. 02*  •  •  • '  "  f  "  f  V%* ^o'  Vo*  •  •  •  ^0^0^ 

00 

X  «xp  [-‘»l«o-‘„?,  »n*l  •  '’o"5o]«o'‘Vo-  •  •  <h 
and  Pg  is  a  double  Fourier  Transform  with  respect  to  both  I  and  C^- 

^2  ^0-^1*  ‘‘2*  -^2 . ^rt>  =j7'^2<^*V^^*  Vo* 

00 

xexp  -UkrV^-i(3o-0l)?o-i2  v"e 

n“  1 

♦‘i:8rtl  ''o5j<15<ls„dv5dv„l„... 
n=l  J 

The  Stochastic  Fourier  Transform  of  Eq.  (A2)  may  now  be  determined  by 
equating  Eq.  (A5)  to  the  sum  of  Eqs.  (A6)  and  (A9).  With  this  and  the  substitution 
ATj  =  •  r"^.  the  following  result  may  be  obtained; 

n^  (2  )^n  /'"/  '^q^^t'^1’^2’--- V*’l^*‘r^l'^2'^2 . ^n'^n^'^^l'^^2- ’  ’  ^^n 


iiTj  .  r.)  < 
t  ‘ 


=  2H^  exp(  +  iki  .  rj)  - 


z  ak2jj  ■  ^‘^y  ak2  -  ^ 


*  ‘  3^'"  liq;  -  ‘  liq-/  <=o> 

^1*  ^2*  ••■ 

X  P2(ki  -  3o*  ^o  "  ^1' *‘2'  *^2*  •  ■  •  ^n*  '^n' ■  ■  •  ‘^^n  ' 


The  first  term  of  the  right-hand  side  of  Eq  (A  10)  has  a  factor  expf  +  ifTi^  •  r^). 
This  implies  that  the  current  may  be  written  in  the  form 

^1’  ^2'  ■  •  •  ^  °  •iq^^l'  ^  exp{  +  ik^  <  r^)  . 

When  this  is  substituted  into  Eq.  (AlO),  the  following  equation  for  may  be  ob¬ 
tained 


S 

(A)  "/■/ V3i*?2'  •••>  Pl<‘‘l-3i.k2-^2*---*^n-^n) 

n-»oo  '  '  L 

-  02‘  ^2‘  '  "  *^^2  *  ■ '  ‘^^n 


2H 


Q  .  Q  3  ,  Q  3 

2  X  akrt..  y  ak^ 


'2x 


y  J 


'Nz'  ^ - ^n^ 


Here,  the  quantity  r2  is  given  by 

^2^®1’  '‘r  ^2'  ^2'  ■  ■  •  ^n^ 


1  --r  ag(Af  .8o)  a  a?(A?‘.8-)a  „  _ 

^  (■i7)2  "l^y  ^ 

X  P2^*^l  *®o*  *®r  ^  exp<ik^  •  Ar^)d 


This  form  of  the  integral  equation  describing  the  current  is  one  that  is  amenable 
to  further  analysis  of  its  component  terms.  The  term-by-term  analysis  leads  to 
determination  of  conditions  for  physical  optics  solutions. 
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